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Abstract 

An exponent (3 which characterises non-equilibrium coarsening processes is 
calculated in a deterministic solvable model of coarsening for a ID q-state 
Potts system. We study how the fraction of sites P which have never changed 
their state, scale with the characteristic domain length {£) . (3 is defined by P 
~ We propose a new model of coarsening that prevents correlations 

from developing between domains thereby ensuring tractability and an exact 
result for any q. 



PACS: 02.50,05.20 
1. Introduction 

Domain coarsening occurs when a system is rapidly quenched from a high 
temperature disordered state to a low temperature ordered state [1]. Do- 
mains of different equilibrium states form on larger and larger length scales 
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thus giving us a system that is forever in nonequihbrium. An important 
observation is that this nonequihbrium coarsening process exhibits dynamic 
scahng at late times, that is, the domain structure at later times is statisti- 
cally similar to those at earlier times except for a global change of scale [1]. 
This implies that at late times the system is characterized by a single length 
scale {i{t)) and all properties of the system can be got in terms of this {i{t)). 
The scaling hypothesis implies that the equal time and two time correlation 
function can be written in scaling form as 

C{r,t) = f{r/{i)) C{r,t,t') = {{i)' / {i))'h{r/ {£)) 

These scaling forms are well supported by experiment. 

Theoretically there have been two main approaches to studying the coars- 
ening process. The first is to study stochastic models such as the Ising model 
with Glauber dynamics and the second is to study deterministic models like 
the Ginzburg-Landau equation for the coarse grained order parameter with 
no noise term. In this case the GL equations are completely deterministic, 
the only randomness being in the initial conditions. This is very different 
from the kinetic Ising models, which are stochastic by definition (even at T 
= 0) and the average of any thermodynamic quantity is not only over ini- 
tial conditions (which is the case for the T = GL equations) but also over 
the histories of evolution. It is generally beheved that in d>2 the stochas- 
tic and deterministic models belong to the same universality class. But in 
ID the situation is quite different. For example for the ID Ising model, 
with Glauber dynamics, at T = 0, evolved from a random initial condition ( 
quenching from high T ) , the domain walls behave as independent random 
walkers who annihilate upon meeting. The average domain size {£) grows as 
t^/^. The equal time and the two time spin-spin correlation functions can 
be exactly calculated [2,3]. On the other hand for deterministic models the 
situation is different. For example in the ID GL model the domain walls in- 
teract via their exponential tails leading to a logarithmic growth law [4] . This 
kind of a difference in time dependence suggests that exponents defined with 
respect to {£) may be more fundamental than those defined with respect to t. 

The correlation function exponents and the dependence of the charac- 
teristic domain length scale on time have been studied extensively [1,2,3,4]. 
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But since this gives no information on the history of evolution of the system, 
the study of "persistence" is important. A simple way of characterizing the 
history of evolution of a system would be to look at, say, the probability of 
a single spin in an Ising system not changing its spin upto a certain time, 
in other words, its "persistence" in clinging to its original state. By study- 
ing the scaling behaviour with (£) of P, the fraction of spins that have not 
changed their state, we can define the exponent (3 by P ^ (i)^^^. For the 
stochastic ID Ising model with Glauber dynamics an exact result for /3 was 
found [5] to be /? = 0.25. The same exponent was studied in the context 
of the ID GL model whose late time dynamics can be mapped, as will be 
explained later, to a deterministic coarsening algorithm which consists of the 
continual removal of the smallest domains. In this case it was found [6] to 
be exactly 0.824... 

The same question can be asked in the context of the q-state Pott's model. 
It turns out that though the autocorrelation exponent and the growth law 
for the characteristic domain length may be the same, f3 varies with q. An 
exact result [5] was recently found for the ID stochastic case, which is an 
exact generalization of the q = 2 case, as 



Noting the difference in value of (3 in the stochastic and deterministic 
case for q = 2, it is only natural to ask how f] varies for general q in the 
deterministic case. For general q, however tractability is lost because, as we 
shall explain later, correlations develop between domains. Nevertheless ap- 
proximate calculations ignoring these correlations have been carried out [7,8]. 

In this paper we present a different model of coarsening in which no cor- 
relations develop between domains. Using this model we can calculate /? 
exactly for any q. For q = 2 our model reduces to the case discussed above. 

2. The Model 

The motivation behind our model is as follows. Let us consider the q 




3 



= 2 Ising case first, in the context of the time dependent Ginzburg-Landau 
equation in ID with no noise ( T = ). It turns out [4] that the late time 
dynamics of this model can be mapped to deterministic equations describing 
the motions of a set of interacting domain walls ( kinks in the wall profile 
). The kinks interact via an exponential attractive interaction. This implies 
that in the limit when the typical separation between domain walls is much 
larger than the width of a wall, the shortest domains collapse instantly as 
compared to the longer ones. This means that the dynamics reduces to a de- 
terministic model where at each time step the shortest domain is found and 
the spins inside it are fiipped. Thus the system coarsens with the continual 
"removal" of the shortest domains. 

Now we can address the question as to what fraction of the spins have 
never fiipped. It is clear that once a domain has been removed the spins 
inside no longer contribute to the fraction of spins that have never fiipped. If 
we call those parts of the line where the spins have fiipped atleast once "wet" , 
and those parts where the spins have never fiipped "dry" , then clearly, when 
a domain is removed, the portion of the line it occupied becomes "wet" . So 
the question about the scaling of the fraction of spins that have never flipped 
can be rephrased as: How does the density of dry regions scale with (£)? We 
define the exponent (3 by saying that the dry region density scales as 
or the dry part per domain goes as 

A straightforward generalization of this procedure to the q - state Pott's 
model can be done. The domains form a random sequence constructed from 
the q available states, no two consecutive domains being the same. Again 
the smallest domains will collapse and can be removed if the states of the 
domains on either side are the same. But if they are different the two walls 
will coalesce to form a new wall at the midpoint. But this routine leads to 
correlations developing between the domains. Instead we consider an ab- 
stract model in which the smallest domain is merged with the domain to its 
right if the states of the domains on either side are different. This ensures 
, as we shall shortly explain, that no correlations develop between domains. 
This leads to tractability and an exact result for any q. 

In this paper we examine how the dry region density scales with the char- 
acteristic domain length scale in the context of this abstract model. 
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The details of the model are as follow. 
At any time we have a sequence of domains constructed from the q available 
states, no two consecutive ones being the same. The smallest domains are 
identified at each time step. If the domains to its left and right are the in 
the same state then it flips to that state ,that is, the three domains merge to 
form a single large domain. Since domains remain uncorrelated, as we shall 
show, the probability of this occurence is On the other hand with a 

probability ( ^^e event that the two domains are in different states ), 
the smallest domain flips to the state of its right side neighbour . In other 
words, it merges with its right hand neighbour leaving the other neighbour 
unaffected. In both cases the part of the line previously occupied by the 
smallest domain becomes wet. 

Now we consider the question of correlations. Suppose we choose N in- 
tervals, the number of distinct arrangements of these on a circle is (N - 1)!. 
Suppose we determistically coarsen this system by picking the smallest do- 
main each time and removing it, till we end up with one interval. Thus (N - 
1)! distinct histories can be created. Now we consider an alternate algorithm 
which consists of picking the smallest interval and then picking two intervals 
at random and combining the three. We iterate this procedure until only 
a single interval remains. This procedure also generates (N - 1)! histories 
which are in one to one correspondence with the histories generated by the 
deterministic algorithm described above. This proves that no correlations de- 
velop between domains during coarsening by merging with domains to both 
the left and right. 

Now consider coarsening the N intervals arranged on a circle by picking 
the smallest interval each time and combining it with the interval to the 
right. This procedure generates (A^ — 1)! histories. As before we consider 
an alternative procedure of picking the smallest interval at each step and 
picking another interval at random and attatching it to the smallest inter- 
val. This procedure also generates (N — 1)! histories again proving that the 
coarsening algorithm does not develop any correlations between the domains. 

So we have proved that the coarsening procedure we use in our model. 
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which is nothing but a combination of the two coarsening algorithms de- 
scribed above, does not allow correlations to develop between intervals. As 
a result we can use the "picking at random" algorithm to compute (3. 

3. Equation for (3 

Our calculations follow the method used by Bray et al [6] for the ID Ising 
model case. We start with random intervals on a hne. Each interval / is 
characterised by its length 1(1) and by the length of its dry part d(I). At each 
time step the smallest domain 7^j„ is picked. As explained before there are 
two possibilities. 

1. Two more intervals /i,/2 are picked at random. The three domains are 
merged to form a single large domain I. This occurs with probability 
The total length and dry parts of / are given by 

l{I) = l{h)+l{Imin) + l{h) (1) 

d{I)^d{h) + d{h) (2) 

2. Another interval Ii is picked at random. The smallest domain is merged 
with it to form a new domain /. This occurs with probability The 
length and dry part of / are given by 

l{I)^l{h)+l{Irmn) (3) 
d{I) = d{h) (4) 

We assume that the lengths of intervals take only integer values and that 
the minimal length in the system is iq. We also assume that there is a very 
large number N of intervals. We denote the number of intervals of length i 
by n,i and also the average length of the dry part of intervals of length i by 
di- Wc denote by primed symbols the values of these quantities after all the 
intervals of length iq have been eliminated, so that the minimal length 
becomes io + 1. The time evolution equations are then given by 

N'^N-2n,,{^)-n:,S^) (5) 
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i.e 



N' ^N-n. 



10^ 



q-r 



(6) 



Similarly we have 



i—2io 



n, = Hi 



J =10 



Tij Tii—j—iQ Q 2 ^n^—iQ\ 



(7) 



H-J-to) 



(8) 



Note that these are vahd for « N which is vahd when ig becomes large. 
Wc assume that after many iterations i.e. when io becomes large, a scaling 
regime is reached, where 



n, 



(9) 



For large io we can treat x = i/io as a continuous variable. Then neglecting 
0(1 /il) we have 



rib- 



and 



^/(^)-^f/(.)-(^)^/(.)-^-.^l (10) 



^iO + 



mr [9{x) - (^) ^g{x) + ^g{x) - X 



io 



io 



(11) 
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Inserting these into the time evolution equations (7) and (8), and using the 
fact that / and g are independent of io for large io we get 



fix) + xnx)+e(x-3)I^ r'f{y)f(x-y-l)dy 



q-lJi 



+ 9{x-2){^)fil)fix-l)=Q 



(12) 



f(l) /■^-2 

(1 - /3)g{x) + xg\x) + 2e{x - 3)^ / g{y)f{x -y- l)dy 



q-lJi 



+ ^(x-2)(i-^)/(lMx-l)=0 
q z 



(13) 



Now we introduce the Laplace tranforms of the functions / and g 

exp{—px)f{x)dx 



exp{—px)g{x)dx 



(14) 
(15) 



Taking Laplace transforms of equations (12) and (13) we get 

p0'(p) = f{l)exp{-p)[{^)^'{p) + {^)Hp) - 1 

pi/j'ip) + pilj{p) = f{l)exp{-p)[2(^——^^ij{p)(f){p) + (^y)^(p) 



(16) 



/(!)- 
(17) 



These are fairly simple first order differential equations in one variable and 
may be solved in a straightforward fashion to obtain the solutions 



1 — exp{—h{p)) 
1 + {^)exp{-h{p)) 



(18) 
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+ ■^(ixp{-h{x))^2exp{h{x)) xl^-^ 



■ 1 + ^ea;p(— -I exp{h{p)) 
where h(x) is given by 



-exp(— (19) 



OO g-t 



(20) 



It may be noted that the constants of integration implied by the form of the 
solutions above, are fixed by the requirement that both and ip go to zero 
as p tends to infinity. 

It remains to fix the constants f(l) and /?. For this we use the following 
expansion 



/ dx = — ln(p) — 7 — 

■^P ^ n=l 



nnl 



(21) 



where 7 is Euler's constant and has a value 7 = 0.577215.... Using this with 
(18) and (20) gives a small p expansion for 



Hp) = 1 - ( 1 



q-l 



pq 



exp 



-/(l)7](l + 0(p)) (22) 



Comparing this with the small p expansion from (14) which is 0(p) = 1 — 
{x)p+ ... we get /(I) = ^ and also the ratio of the mean domain length to 

the minimum length as (x) — (^zj-)^^- 



by 



The exponent /3 can be determined in a similar fashion. We define r(p) 



r{p) = h(p) + ln(p) = -7 - ^ 



n=l 



nnl 



(23) 



Using this we can rewrite (19) as 

■1 + ^exp{-r{x)).2exp{r{x))x^ 



Jp 



1 + -^exp{—r{p)) J exp{r{p)) pf^ 



—exp{—x)dx (24) 
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From this we may get the small p form of ip{p) as 



f^ + f^-^V-^Bfe/J) (25) 



where 

roc ^ ^ d 



roo c\ r / r(x) —r(x) \ 2l 

B(p,l3)= x^-'-^\e-^(e^ + xe^) ]dx (26) 

Jp dx L V ? — 1 ^ 

It may be easily shown that 

B{p, p) = 5(0, p) + 0{p^+^) + Oip'+I^) + ... (27) 

Thus (25) reduces to 

^(l + S(0, P)p'-^ + 0{p) + 0{/) + ...) (28) 

We now compare this expansion to the direct expansion of ■0(p) obtained 
from (15), namely 'ip{p) = dxg{x){l — p{x) + 0(p^)). If we require the 
function g{x) to have a finite first moment then we have to have B{0, (5) — 
i.e 

^i^-^—L-(^e^ + -±-^xe^) ]dx = Q (29) 



" dx 



This condition determines /3 for us. For general q the /3 can be determined 
by numerical methods from the above condition (29). We have computed /? 
for sample values of q (see table). /3 for q = 2 agrees with Bray's value [6] 
upto the 5 decimal places we have computed. This is what one would expect 
because, as we mentioned before, our model reduces to Bray's model for q = 
2. We also notice a monotonic decrease in /3 as q increases. This is also what 
one would physically expect if we look at the dry part per domain which 
scales as (i)^. As we go to higher values of q, coarsening proceeds almost 
exclusively by the smallest domain merging with one nearest neighbour. This 
implies that, in the scaling regime, the dry part per domain remains almost 
constant leading to lower and lower values of (5, which indeed, is what we 
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find. It can also be shown analytically by considering the small p form of 
ip{p) for the q— > oo case that P = 0. 

4. Conclusion 

Wc have studied a deterministic model of coarsening for the zero temper- 
ature dynamics of a q-state Pott's system in ID following a rapid quench. 
This model does not allow correlations to develop between domains thus ren- 
dering it exactly solvable. We have determined the persistence exponent /5 
exactly for all q. For q = 2 we obtain f3 = 0.824... in complete agreement 
with [6]. As q increases we find that /3 decays monotonically to zero. 
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Table 1: /3 for various q 



q 


/9 


2 


0.82492... 


3 


0.68092... 


4 


0.58178... 


5 


0.50940... 


7 


0.40999... 


10 


0.31905... 


50 


0.08321... 
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